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ABSTRACT 

The oblique geometry of pulsar wind termination shock ensures that the Doppler 
beaming has a strong impact on the shock emission. We illustrate this using the recent 
relativistic MHD simulations of the Crab Nebula and analysis of oblique shocks. We 
also show that the observed size, shape, and distance from the Crab pulsar of the 
Crab Nebula inner knot are consistent with its interpretation as a Doppler-boosted 
emission from the termination shock. If the electrons responsible for the synchrotron 
gamma-rays are accelerated only at the termination shock then their short life-time 
ensures that these gamma-rays originate close to the shock and are also strongly 
effected by the Doppler beaming. As the result, bulk of the observed synchrotron 
gamma-rays of the Crab Nebula around 100 MeV may come from its inner knot. This 
hypothesis is consistent with the observed optical flux of the inner knot, provided 
its optical-gamma spectral index is the same as the injection spectral index found 
in the Kennel & Coroniti model of the nebula spectrum. The observed variability of 
synchrotron gamma-ray emission on the time scale of wisp production can be caused 
by the instability of the termination shock discovered in recent numerical simulations. 
Given the small size of the knot, it is also possible that the September 2010 gamma- 
ray flare of the Crab Nebula also came from the knot, though the actual mechanism 
remains unclear. The model predicts correlation of the temporal variability of the 
synchrotron gamma-ray flux in the Fermi and AGILE windows with the variability of 
the un-pulsed optical flux from within 1" of the Crab pulsar. 
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1 INTRODUCTION 

The Crab Nebula has been a source of intriguing discov- 
eries and served as a test bed of astrophysics for decades. 
This is one of the best studied objects beyond the Solar sys- 
tem. It has been observed at all wavelengths, from radio to 
very high energy gamma-rays. Its non-thermal emission be- 
low £p h ~ 500 MeV is a synchrotron emission of relativistic 
electrons in the nebula magnetic field and above £S h it is 
the inverse Compton emission of the same electrons. The 
emitting electrons are accelerated up to PeV energies, indi- 
cating that the acceleration mechanism is very efficient. The 
source of energy is the ultra-relativistic magnetic wind from 
the pulsar ( |Rees fc Gunn||1974| |Kennel fc Coroniti||1984a[ ), 
but the actual mechanism of particle acceleration is still a 
mystery. The main candidates are the diffusive shock ac- 
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celeration at the wind termination shock, the second-order 
Fermi acceleration in the turbulent plasma of the nebula, 
including secondary shocks, and the magnetic reconnection 
events. 

Compared to the highly filamentary thermal emission, 
the non-thermal emission is relatively featureless. Yet, it was 
discovered already in 1920 that fine and dynamic "wisps" 



are somehow produced in the center of the nebula (Lamp- 
land|1921| [Scargle|1969"| . Later, the X-ray observations dis- 



covered the famous jet-torus structure in the inner nebula 
(Weisskopf et al.|2000"}|Hester et al.|2002|), and the high reso- 



lution optical observations with Hubble Space Telescope re 
vealed fine sub-arcsecond structure of the non-thermal emis- 
sion, including few optical knots ( Hester et al.|1995 20021. 

The synchrotron life-time of electrons emitting in op- 
tics is comparable to the dynamical time-scale of the nebula, 
and this makes it difficult to spot the exact locations of the 
particle acceleration cites. In gamma-rays below £p h , where 
the life-time becomes short, the angular resolution of the 



2 Komissarov & Lyutikov 



telescopes is insufficient to see the nebula structure. How- 
ever, the observations indicated variability of the gamma-ray 
emission in the 1-150 MeV range ( Much et al.|1995||de Jager| 
et al. 1996 ) on the time scale around one year, de Jager et al.l 



I 1996 ) proposed that this emission could originate from the 



variable optical features seen with HST in the polar regions 
of the inner nebula, in particular the so-called "anvil" . Vari- 
ability on a similar time-scale has been recently discovered 
in the X-ray emission ( |Wilson-Hodge et al.||2010| ). 

In September 2010, AGILE collaboration reported a 
three- fold increase of the gamma-ray flux (>100 MeV) from 
the direction of the Crab Nebula ( jTavani et al.|20l"Tj ) , which 
was immediately confirmed by Fermi LAT collaboration, 



who reported a six-fold increase of the flux (Abdo et al. 
2011). The flare continued for four days, September 18-22, 
after which the gamma-ray flux returned to the pre-flare 
level. Fermi also reported that the pulsed emission of the 
Crab pulsar remained unchanged during the flare, suggest- 
ing that the flare originated in the Nebula. Jodrell Bank ra- 
dio timing observations of the Crab pulsar showed no glitch 
during the flare, supporting this conclusion (ATel#: 2889). 
At the same time, INTEGRAL reported no detection of the 
flare during Sep 19 in the 20-400 keV window (ATel#: 2856) 
and Swift/BAT did not see any significant variability dur- 
ing the gamma-ray flare in the 14-150 keV range (ATel#: 
2893). Swift also reported no evidence for active AGN near 
the Crab, suggesting that the Crab itself is responsible for 
the flare (ATel#: 2868). ARGO-YBJ collaboration reported 
a significant enhancement of the very high energy emission, 
around 1 TeV, from the Crab nebula during the AGILE- 
Fermi flare (ATel#: 2921). However, this has not been con- 
firmed by VERITAS and MAGIC collaborations (ATel#: 
2967,2968). This discovery seems to have given credit to 
another event, detected in February 2009, which lasted for 
approximately 14 days, during which the gamma-ray flux 



increased by a factor of three or four (Tavani et al. 2011 
Abdo et al.|2011 " |. On the SED plots the flares appear as an 



extension of the synchrotron component further out towards 
the higher energies, up to 1 GeV for the September flare and 
a bit less dramatic for the February flare. 

The short duration of these flares suggests that their 
source is rather compact. Unfortunately, no high angular 
resolution observations of the nebula were carried out dur- 
ing the flares. The Crab Nebula images from Chandra and 
HST, obtained after the September flare, have not revealed 
anything especially unusual (ATel#: 2882, ATel#: 2903). 
They do show a change in the structure of the nebula wisps 
compared to previous observations, carried out years ear- 
lier. However, the large length scale of these wisps shows 
that they can hardly be a source of the flares. The Chan- 
dra images also show a significant change in the position 
of one of the jet knots, which apparently had moved about 
3" towards the pulsar. This may be more significant as this 
feature is more compact. 

Even more potent source of the flares could be the mys- 
terious "inner knot", discovered only 0'.'65 away from the 
Crab pulsar along the jet direction. This knot, named as 
"knot 1" in |Hester et al.| ( |1995[ ), is the brightest and most 
compact persistent feature of the Crab Nebula. It is seen 
at more or less the same location in many observations, 
both with space and ground based telescopes with adap- 
tive optics, which followed its discovery. It is extended and 



elongated, with the main axis perpendicular to the jet. Its 
length and width are ipx ~ 0'.'50 and ipn ~ 0'.'16 respectively 
( jHester et al.| [l995). A number of recent optical and infra- 
red observations of the inner Crab Nebula focused on the 



knot variability. Melatos et al. ( 2005 1 reported no variabil- 



ity on the timescales from 1 kilosecond to 48hr. Tzi amtzis| 
et al. ( 2009 1 compared the measurements separated by two 



and half months and found no significant difference too. On 



the other hand, Sollcrman (2003) analyzed the HST archive 
data and discovered flux variations on the level of 50% over 
the period of 6 year^] They also reported possible random 
displacements of the knot on the scale of O'.'l. More recently, 
Sandberg fc Sollerman| ( |2009 1 reported twice as higher flux 
from the inner knot in 2003 compared to the earlier observa- 
tions in 2000. Thus, we tentatively conclude that the inner 
knot does show significant variations of its flux, and pos- 
sibly location, on the time-scale comparable to that of the 



gamma-rays variability reported by de Jager et al. ( 1996 1 



These data suggest to consider the possibility that the inner 
knot may be a strong source of gamma-rays, both during 
and between the flares. 

The synthetic maps of synchrotron emission from the 
Crab Nebula, based on numerical relativistic MHD sim- 
ulations, reveal a bright compact feature strikingly simi- 
lar in its appearance and location to the inner knot (see 
Figure [2j|. In Komissarov & Lyubarsky (20041 this feature 
was identified with the Doppler-boosted emission from the 
high-velocity flow downstream of the oblique termination 
shock of the pulsar wind. The more recent simulations of 
the Crab Nebula, which had a significantly higher resolu- 
tion, discovered strong variability of the termination shock, 



involving dramatic changes of its shape and inclination ( Ca- 
mus et aL]|2009 1. This discovery suggests that the gamma- 



ray variability may be related to the changes in the Doppler 
beaming of the post-shock emission, associated with this 
structural variability of the termination shock. |Volpi et al.| 
( 2008 ) have already addressed the issue of variability of both 
the synchrotron and inverse Compton emission of the Crab 
Nebula in their computer simulations. They reported the 
strongest variability in the part of the gamma-ray window 
where the synchrotron emission is still dominant over the 
inverse Compton emission, and on the time-scale compara- 
ble to that found by ( |de Jager et al.|[l996[ ). However, they 
could not identify the source of variability with any partic- 
ular feature in their numerical solutions. 

Here we present new arguments in favour of the inter- 
pretation of the inner knot as a Doppler-boosted shock emis- 
sion, and investigate whether the inner knot can be a strong 
source of variable gamma-ray emission. In Section[2]we anal- 
yse the observed shape and location of the inner knot and 
show that they are consistent with this interpretation. In 
Section [3] we argue that at least a significant fraction of the 
synchrotron gamma-ray emission of the Crab Nebula orig- 
inates from the inner knot. The key factors are the short 
cooling time of electrons and the strong Doppler beaming 
of the emission originated in the vicinity of the termina- 
tion shock. In Section [4] we discuss the possible connection 
between the observed variability of the gamma-ray emission 



1 For some reason the actual observational data have not been 
presented in this paper. 
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Figure 1. Interpretation of the inner knot as the Doppler-boosted 
emission from the high velocity flow located downstream of the 
oblique termination shock. The termination shock is shown by 
the solid line. The dashed line is the line of sight. 



the proper magnetic field in this region is strongest. Several 
factors are responsible for this result. Although the the mag- 
netic field of dissipationless wind behaves as B$ ex aiaO/r 
l Michel|197 3 Bogovalov 1999), in the simulations this func- 
tion was multiplied by 1 1 — 29/tt\ , in order to account for the 
magnetic field dissipation in the striped wind zone. Thus, the 
wind magnetic field peaked at 9 ~ 57° instead of 9 ~ 90°. 
More important, however, is the axial compression of the 
nebula by the magnetic hoop stress, which leads to the total 
pressure downstream of the termination shock to be signif- 
icantly higher at small polar angles. As the result, the arch 
shock is pushed closer to the pulsar, leading to a stronger 
upstream and hence downstream magnetic field. Second, the 
emission from the upper arch shock is Doppler boosted. In- 
deed, as one can see in the left panel of Figure [2] the Lorentz 
factor of the flow downstream of the arch shocks is quite 
high. 

The analysis of oblique relativistic MHD shocks given 
in Appendix [X] shows that in the case of ultra-relativistic 
cold upstream flow the downstream Torentz factor is 



and the non-stationary shock dynamics discovered in numer- 
ical simulations, in particular the role of the variable Doppler 
beam orientation. Our conclusions are given in Section [5] In 
Appendix|A]we present the analysis of oblique MHD shocks. 
Its results allow us to determine how high the Lorentz factor 
downstream of the Crab's termination shock can be and to 
confirm the results of numerical simulations. 



2 ORIGIN OF THE INNER KNOT 

The first numerical simulations of the Crab Nebula by 



72 



( Komissarov fc Lyubarsky 2003 2004 Del Zanna et al.|2004 
|Bogovalov et al.||2005| ) had rather low numerical resolution 
and imposed rcflcctional symmetry in the equatorial plane. 
In these simulations the termination shock appeared as a 
more or less stationary complex structure, which in fact in- 
cluded several different shocks. These components, named as 
arch shocks, rim shocks, and the Mach belt in Komiss arov fc| 
|Lyubarsky|([2004l), are shown in Figure [I] The latest simula- 
tions by Camus et al. ( 2009 1 had much higher resolution and 



the computational domain included the whole range of the 
polar angle, 9 £ (0, n). In these simulations, the structure of 
termination shock appeared highly distorted and dynamic, 
but as one can see in Figure [2j these individual components 
were still identifiable. 

The right panel of Figure [2] shows the typical distri- 
bution of the synchrotron emissivity as measured in the 
pulsar/observer frame, found in the latest numerical sim- 
ulations]^] One can see that it is strongly enhanced in the 
vicinity of the upper arch shock of the termination shock 
complex. There are two reasons for this enhancement. First, 



2 The details of these numerical simulations and the radiation 
transfer computations are described in |Camus et al.| | |2009| |. Here 
we only note that the radiation model assumes that relativistic 
electrons with the power law energy spectrum N(£ e ) oc £ 22 , 
terminated at £™ ax = 1 PeV, are injected at the termination 
shock, to be more precise at the arch shock and the Mach belt, 
and then they evolve subject to advection and synchrotron and 
adiabatic energy losses. 



(1) 



where <5i is the angle between the upstream velocity and the 
shock plane and 



X = 



1 + 2<7i + ^16cr 2 + 16cti + 1 

6(l + <Tl) ' 



(2) 



where <ti = B 2 /4npic 2 is the magnetization parameter of 
the upstream flow. For o\ ^> 1 this yields 
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sin Si 



and for 
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2^2 



(1- 



1 ^ 1 
2 ai) sinV 



(3) 



(4) 



In the simulations ai varies with the polar angle between 
and 0.05, with the volume averaged value (cri) ~ 0.014, 
and thus the latter limit applies. On can see that for small 
shock inclination angles the Lorentz factor can indeed be 
quite high. Even higher values are expected for high-sigma 
pulsar wind. 

In addition to having high Lorentz factor in the down- 
stream flow, the upper arch shock is inclined at the angle 
of ~ 60° to the polar axis, near the axis. Observations of 
the inner Crab Nebula suggest that the angle between the 
line of sight and the symmetry axis of the nebula is also 
close to 60° ( jWeisskopf et al.||2000] ). Thus, the upper arch 
shock is well aligned with the line of sight, resulting in strong 
Doppler-boosting of its emission. This is schematically illus- 
trated in Figure [l] The left panel of Figure [3] shows the 
synthetic optical synchrotron image of the inner part of the 
simulated PWN at the time corresponding to the age of the 
Crab Nebula. One can see prominent wisps and a bright 
knot located very close to the origin, where the projected 
image of the pulsar would appear if it was included in the 
emission model. In the simulations, there is no emission from 
the pulsar wind as the emitting electrons are injected at the 
termination shock only. Thus, all the fine features of the 
synthetic synchrotron images, including the inner knot, are 
produces inside the nebula. 
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Figure 2. Termination shock in numerical simulations of |Camus et al,| ( |2009[ ) . The left panel shows the flow Lorentz factor. In this plot 
the wind zone has been cut off (The thin line at r ~ 0.7 is an artifact of this procedure, which conveniently indicates the location of the 
Mach belt). The right panel shows the observed synchrotron emissivity in the optical range, log 10 J„ l0 b, in the plane which includes the 
line of sight and the symmetry axis. The angle between the line of sight and the symmetry axis is 60°. 
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Figure 3. Synthetic images of the inner Crab nebula in optics. The panel shows the proper synchrotron image, log 10 I v , where I v is 
the intensity of radiation. The angle between the line of sight and the symmetry axis is 60° . The right panel shows how the image would 
look if the the Doppler beaming was not taken into account. 



It is interesting that both in the synthetic and the real 
optical images of the Crab Nebula the termination shock 
is not clearly identifiabl^] The inner wisps give away its 

3 The Chandra X-ray image of the inner Crab Nebular is overall 
similar to the HST image l |Hester et aLl|2002| . There are how- 
ever some noticeable differences. In particular, the X-ray image is 
much knottier and some of these X-ray knots arrange in a sort of 



size but that is about it. Several factors contribute to this 
effect. Firstly, the inner cavity filled with the pulsar wind 

ring, called the "inner ring" , around the pulsar. This ring is often 
identified with the termination shock. The fact that the ring is 
relatively symmetric, in contrast with the optical image, is hard 
to explain in our model and indicates that some other factors 
have to be included. 
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Figure 4. Basic geometric parameters of the inner knot. 



is small compared to the size of the nebula, making the 
central brightness reduction rather weak. Secondly, the re- 
gion of enhanced proper emissivity around the termination 
shock forms a geometrically thin distorted shell. Because of 
this, the observed emission is strongly enhanced in places 
where the line of sight is tangent to the shell surface. This 
leads to the appearance of several bright rings on synthetic 
maps where the Doppler beaming is not included (see the 
right panel of Figure [3J. The Doppler beaming leads to in- 
creased emissivity in the part of such a ring where its plasma 
flows towards the observer, and decreased emissivity where 
it flows away, turning the ring into an arc (see the left panel 
of Figure [3} . Because of the non-spherical shape of the ter- 
mination shock, the jet base is located much closer to the 
pulsar compared to the Mach belt radius, making the wrong 
impression that the jet originates from the inside of the ter- 
mination shock and suggesting that it is produced by the 
pulsar. 

Comparing the images presented in Figure [3] one can 
see how some features brighten up and others get dimmer 
because of the Doppler beaming. In particularly, the bright 
knot in the center of the left image is certainly Doppler- 
boosted. |Komissarov fc Lyuba rsky (2004) proposed that this 
synthetic knot is a counterpart of the Crab's inner knot. In- 
deed, like the inner knot it is positioned on the jet-side of 
the nebula at the base of the jet (only in projection) and it 
is elongated in the direction perpendicular to the jet ( Hester 



|et al.|1 995 ) . Other synthetic jet knots, which are seen in Fig- 
urc[3] have more or less the same brightness in both these im- 
ages, indicating that the Doppler beaming is not that impor- 
tant. They indeed originate at the base of the polar jets. The 
jet Lorentz factor is variable but on average it is rather low, 
7j ~ 1.5. Combined with the large viewing angle, 60°, this 
explain why the Doppler effect is rather weak. Phenomeno- 
logically, these jet knots are created via unsteady inhomo- 
geneous axial magnetic pinch, which is responsible for the 
jet formation. The flow towards the polar axis, which feeds 
the jet, is highly inhomogeneous with patches of strong and 
weak magnetic field, resulting in strong spatial variations of 
the magnetic hoop stress and highly variable jet dynamics. 
The non-linear sausage-mode of the magnetic pinch insta- 
bility could be another contributor to the jet variability. 

The fact that the inner knot of the Crab Nebula is so 



compact and that it appears at the base of the Crab jet 
are more suggestive of some jet feature rather than the ex- 
tended termination shock. However, the knot shape and size 
are nicely explained in the shock model. Denote as tp p the 
observed angular distance of the knot from the pulsar, as d n 
the linear distance to the nebula, and as the linear dis- 
tance between the pulsar and the point on the shock surface 
where the line of sight is tangent to the shock. The angle 
between the shock surface and the upstream velocity vector 
denote as <5i (see Fig|4|. Then in the small angle approxi- 
mation 



Si = ip p (d n /d k ) 



(5) 



Using EqjT] we can now find the Lorentz factor of the post- 
shock flow 



72 = 



k 

5i 



(6) 



4>D 



(7) 



and the beaming angle 

_L = §L 

72 U ' 

where f a = (1 — x(o"i)) _1//2 - This allows us to find the trans- 
verse angular size of the knot 



dk<f>D 
d n 



Jo 



For <7i < 1, one has f a ~ 1 and, thus, 

tp-L - Ipp, 



(8) 



(9) 



which is in excellent agreement with the observations of the 
inner knot. The same argument shows that HST knot 2, 
which has similar size but is located much further out from 
the pulsar, cannot be the Doppler-boosted part of the termi- 
nation shock, unless the pulsar wind magnetization o\ ^> 1. 

The other size of the inner knot, ip\\, is determined by 
the thickness of the post-shock plasma flow and the shock lo- 
cal curvature. The flow thickness at the distance du from the 
pulsar can be estimated as &2dk, where 82 is the downstream 
angle between the velocity vector and the shock plane. Since 
the distance of the shock from the pulsar across the line of 
sight is Sidk, this yields ip\\ = (<$2/<5i)Vv Using Eq A19 to 
evaluate 82/81 and Eq[9]we then find that if 1/)^ is fully de- 
termined by the thickness of the post-shock flow then 
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(10) 



the termination shock. A number of different ideas have been 



in excellent agreement with the observations again. 

The shock curvature would lead to a finite linear width 
of the knot in projection on the plane of the sky even if 
the postshock flow was infinitely thin. This width can be 
estimated as R C {1 — cos(<^d/2)), where R c is the local shock 
curvature radius. For the small angles this yields 



(Rc/d k ) 

872 u 



(11) 



Thus, unless {Rc/dk) > 872/0-, we have ip\\ < ip±, in agree- 
ment with the observations. Numerical simulations show 
that normally Rc/dk < 10. 

The transverse angular size of the knot could have been 
used to infer the Lorentz factor of the post-shock flow if 
we knew the distance d k between the knot and the pulsar. 
Indeed, combining Eqs. [6] and [5] one finds that 

du J_ 

dn 1p± ' 



72 



Since we only know that dk cannot exceed the radius of the 
termination shock, this equation only allows us to find the 
upper limit 



72 < 



20, 



where i/i ts — 10" is the angular size of termination shock 
inferred from the observations ( Hester et al.|2002 l. 

If the inner knot is indeed a part of the termination 
shock then its spectrum can be used to infer the proper- 
ties of the particle acceleration at the shock. Several groups 
have carried out optical and near-infrared observations of 
the knot in recent years. Unfortunately, their results do not 



quite agree. Sandberg fc Sollerman| ( |2009[ ) reported the opti- 
cal spectral index a = 1.3, assuming that I v oc v~ a . On the 
other hand, the data presented in their Fig. 2 suggest a much 
flatter near-infrared spectrum, with a ~ 0.3. It is difficult to 
see how the synchrotron mechanism can accommodate such 
a large variation of spectral index within only one decade 
of frequency. Perhaps, these measurements suffer from large 
systematic errors. The proximity of the knot to the pulsar 
could be one of the complications. According to the data 



obtained by Melatos et al. ( 2005 1 the near- infrared spectral 



index of the inner knot is a = 0.78. Unfortunately, the ac- 
curacy of this measurement is not given, whereas for other 
features the error is given and it is about ±0.13. The proxim- 
ity to the pulsar suggests that for the inner knot the error is 
higher. Finally, |Tziamtzis eTaL] p009]) give a = 0.63 ± 0.02 



for the optical emission of the inner knot. The very small 
error indicates that this is the most accurate measurement 
to date. This result is in excellent agreement with the value 
of the injection spectral index inferred by Kennel & Coroniti 
( |1984[ > via model fitting of the integral spectrum of the Crab 
Nebula. 

Although, the MHD model in general, and the recent 
numerical MHD simulations in particular, have enjoyed a 
lot of success in explaining the properties of the Crab Neb- 
ula, as well as other PWN, it is by no means problems free. 
The so-called cr-problem is its main weakness. It is not clear 
how exactly the pulsar wind turns from being Poynting- 
dominated near the pulsar to kinetic-energy-dominated near 



put forward but the issue is far from settled (Lyubarsky 



Kirk|200l| |Kirk fc Skjaeraasen||2003l |Lyubarsky||2003b|al 
Arons 2008). An alternative model, where the flow remains 



Poynting-dominated even inside the nebula, has been put 



forward recently (Lyutikov 20101. As the result, it is not 



clear as to what model of the pulsar wind should be used in 
setting the inflow boundary condition in the MHD simula- 
tions of PWN. Moreover, so far only two-dimensional sim- 
ulation have been carried out, which leaves unexplored the 
effects of non-axisymmetric instabilities on the nebula struc- 
ture and dynamics ( |Begelman| 1998[ ) . 



How and where the emitting particles are accelerated 
is also debated. In the simulations, it was assumed that the 
synchrotron electrons are accelerated only at the arch shock 
and the Mach belt, but the acceleration can also occur at 
the rim shocks and in the turbulent interior of the nebula. 



Given these uncertainties, it is not surprising that there 
are still some significant quantitative differences between the 
theory and the observations. For example, the observations 
reveal comparable isotropic luminosities of the inner knot 
and the brightest wisps (Hest er et al.fl9 95). In contrast, in 
the synthetic synchrotron images the knot is several times 
brighter. This is illustrated in Figure [2] which shows the 
intensity of radiation using logarithmic scaling. One possi- 
ble reason for this is the excessive axial compression of the 
nebula via the hoop stress of the azimuthal magnetic field, 
caused by the condition of axisymmetry. In fact, the syn- 
thetic optical images of the Crab Nebula show strong global 
enhancement of the surface brightness along the symmetry 
axis, which is not present in the images of the Crab Neb- 
ula. The observed large scale kink of the Crab jet indicates 
that some kind of kink-mode instability significantly reduces 
the degree of symmetry in the polar regions of the nebula. 
This is important as, in addition to the anisotropic power 
distribution in the pulsar wind, the overall geometry of the 
termination shock is also influenced by the pressure distri- 
bution inside the nebula. The enhanced pressure near the 
axis pushes the cusp of the arch shock further down towards 
the pulsar, making the termination shock less spherical com- 
pared to what it would be in the case of uniform pressure 
distribution inside the nebula. This leads to higher magnetic 
field and number density of emitting particles, and hence in- 
creased volume emissivity near the shock cusp. In addition, 
the shock geometry determines the velocity distribution and 
hence the Doppler beaming. 



In three dimensional simulations of the Crab Nebula, 
the almost perfect alignment of the arch shock with the line 
of sigh may no longer be found. Less squashed along the po- 
lar direction, the termination shock would have lower proper 
emissivity near the cusp region. The Doppler-boosting of the 
knot emission will be reduced as well, and not only because 
of the less perfect alignment of its Doppler beam with the 
line of sight. More spherical shape of the termination shock 
would also lead to higher shock inclination angles and lower 
Lorentz factor of the downstream flow. 
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Figure 5. Synthetic synchrotron map of the Crab Nebula at 
100 MeV. The image shows log 10 /„ in arbitrary units. 



3 WHY THE INNER KNOT CAN BE A 
STRONG SOURCE OF GAMMA-RAY 
EMISSION 

Suppose that the termination shock is the main accelera- 
tion site of gamma-ray emitting electrons. The synchrotron 
cooling timescale is 



tc. 



■" D ' n (lis) 



-3/2 



( £ph,ob 

V 100 MeV 



-1/2 



days, 



where £ p h,ob = D£ P h is the observed energy of photons emit- 
ted at the energy £ p h in the fluid frame of the downstream 
plasma, 



D = 



7(1-/9 cos ( 



is the Doppler factor, B is the magnetic field strength as 
measured in the fluid frame, and <j> is the angle between 
the line of sight and the velocity vector of the plasma 
bulk motion. Fitting of the Crab Nebula spectrum with the 
synchro- Compton model yields the typical B ~ 100 — 200 /j,G 
( Horns fc Aharonian||2004| [Abdo et al.||2010[ ). However, in- 
dividual bright features can have stronger magnetic field. 
In particular, |Hester et al.| ( |1995[ ) give the equipartition 
B ~ 2.5 x hHg for the inner knot. Moreover, the numeri- 
cal simulations show that near the arch shock the magnetic 
field can be significantly higher, up to ten times, than the 
volume averagecQ Thus, even for the Doppler factor as high 
as D = 10, the cooling length scale of electrons emitting at 
£ P h ~ 100 MeV is likely to be small compared to the ter- 
mination shock radius, which is about ~ 10" or ~ 120 light 



days ( Hester et al.| 2002) in linear scale. Thus, the gamma- 
ray emitting region must be located very close to the shock. 



4 This is in contrast with the one-dimensional MHD model where 
the magnetic field increases with the distance from the termina- 
tion shock (Ke nnel &: Coroniti|1984a l 



Downstream of the Mach belt the magnetic field and 
the number density of the emitting particles is significantly 
lower compared to the arch shock, and for this reason the 
synchrotron emissivity is low as well (see Figure |2j. Thus, 
it is the wind plasma which has just passed through the 
arch shock which is likely to be the main emitter of the 
synchrotron gamma-rays. This plasma flows with relativis- 
tic speed and is subject to strong Doppler beaming (see Fig- 
ure]^]). This results in boosting of the emission from the part 
of the flow where the velocity vector is close to line of sight, 
the inner knot region, and dimming of the emission from 
other parts, where the viewing angle exceeds 2/72. Thus, a 
significant fraction, if not most, of the observed synchrotron 
gamma-ray emission of the Crab Nebula may originate from 
its inner knot. In fact, quick inspection of our numerical so- 
lutions shows that at 100 MeV the inner knot is essentially 
the only feature in the sky (see Figure[5|. However, given the 
uncertainties of the numerical model, one cannot exclude a 
contribution from few brightest wisps. 

In order to test this idea against the observations one 
can compare the observed flux from the Crab's inner knot 
in optics with the observed flux from the whole of the Crab 
Nebula at 100 MeV. Given the small light crossing time of 
the knot compared to the synchrotron cooling time even at 
100 MeV, its synchrotron electrons must still have the en- 
ergy spectrum which is very close to the one produced by 
the shock acceleration mechanism. At v = 3.76 x 10 14 Hz 
the de-reddened flux fr om the inner knot is F v ~ 1.6 x 
1 (jTziamtzis et al|2009|). At 100 MeV 



10"^ergs _i cm" J Hz 
the observed flux is F, 



1.7 x 10 _dj! ergs~ i cm~' ! Hz" 



(|Abdo et al. 20101. The corresponding two point spectral 



index is a ~ 0.64. This is indeed the value of the injection 



spectral index inferred by Kennel & Coroniti (19841, and 



measured in optics by Tziamtzis et al. ( 2009 \ 



There exists an upper limit on the energy of synchrotron 
photons, which is independent on the details of the accelera- 
tion mechanism ( Lyutikov|2010 1. If the accelerating electric 
field E is a fraction rj < 1 of the magnetic field then the rate 
of energy gain can be estimated as 



= eEc — rjeBc. 



(12) 



acceleration time scale r acc 
-1 



d£ c 
dt 

The corresponding 

£c/(d£ c /dt) = (t^ojb)^ 1 , where cjb = ceB/£ a is the rela- 
tivistic Larmor frequency. The energy loss rate due to syn- 
chrotron emission 
d£ a 
dt 

where c 2 = 4e 4 /9m 4 c 7 and we also assumed effective pitch 
angle scattering, grows with the electron energy. The balance 
of energy gains and losses yields the maximum energy, which 
can be reached by the accelerated electrons 

= (riec/c 2 B) 1/2 . 



= -c 2 B 2 £c, 



(13) 



(14) 



The characteristic energy of the synchrotron photons emit- 
ted by the electron of energy £ , ™ ax in the magnetic field of 
strength B, 

£ P T X = Cl B(£ri 2 = ~r) ^ = 236 V MeV , (15) 
167T e z 

where c\ = 3e/i/47rm 3 c 5 , does not depend on the magnetic 
field strength. This is the utmost upper limit, which may be 
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impossible to reach in practice. For example, dc Jagcr ct al. 
( 1996 1 give an almost ten times smaller value for £^ x for 
the shock acceleration. 

In fact, the synchrotron component of the Crab Nebula 
spectrum becomes very steep above lOMeV, and can be 
fitted with the function 



F v oc v a exp(— hiy/£p h ), 



(16) 



where the cut-off energy ££ h ~ 100 MeV |Abdo et al.|2010 | 



This is surprisingly close to our value of fcpif^ However, the 
above limit applies only in the frame of emitting plasma. If 
the plasma is moving with relativistic bulk speed relative 
to the observer then it has to be multiplied by the Doppler 
factor in order to obtained the corresponding observed pho- 
ton energy. For 7 > 1 the maximum value of the Doppler 
factor is D m ax — 27, and thus even for the rather moderate 
postshock value of 72 — 5 the synchrotron cutoff energy can 
be increased by a factor of ten. This has been used to argue 
that the observed synchrotron emission of the Crab Nebula 



with £ p h > 100 MeV originates in relativistic flow ( Lyutikov 



20101. Here, this argument can be refined to support the 



inner knot as the source of this gamma-ray emission be- 
cause we are now almost certain that its emission is indeed 
Doppler-boosted. 

It seems reasonable to expect the gyroradius of electrons 
accelerated at the termination shock to be below the radius 
of the termination shock. An effective Doppler boosting may 
also require the gyroradius to be below the transverse size of 
the fast flow downstream of the arch-shock. Only in this case 
one can firmly conclude that the electrons are advected with 
the flow. For the electrons emitting synchrotron photons of 
the energy £ p h in the comoving frame, the gyroradius radius 
is 



£ph 

cie 2 B 3 



1/2 



(17) 



Our numerical simulations show that the magnetic field mea- 
sured just downstream of the arch shock is significantly 
higher than the volume averaged one, up to about one order 
of magnitude. This indicates that B = 10 -3 G may well be 
typical for this region. Then the typical gyroradius radius of 
the electrons is 



1.6 / £ph,ob 

D 1 /' 2 V 100 MeV 



1/2 



B 
10" 3 G 



-3/2 



light days, (18) 



which is significantly less than the termination shock radius 
and even below the size hi ~ 2 light days of the inner knot, 
which is an observational indicator of the thickness of the 



fast post-shock flow (see the discussion leading to Eq 10 1 



Given this result, we conclude that the observed emission 
up to £ p h,ob = 1 GeV can be the Doppler-boosted emission 
produced by the electrons accelerated at the termination 
shock. 



4 NATURE OF THE GAMMA-RAY 
VARIABILITY 

The strong variability of the termination shock discovered in 
the numerical simulations by | Camus et~aT1 ( [2009l ) and associ- 
ated with the wisp production could be behind the observed 



variability of the gamma-ray emission from the Crab Neb- 
ula. Figure [6] shows the intra-year variability of the inner 
knot at £ p h = 100 MeV based on the results of these sim- 
ulations. The epoch corresponds to the present age of the 
Crab Nebula. Within this particular period the total flux 
changes were limited by ~ 27%. (Unfortunately, most of the 
simulation data is now lost and we cannot comment on the 
statistical significance of this result.) The data also indicate 
noticeable changes in the knot appearance and small changes 
in its location. 

The mechanism of this shock variability is not very 
clear. It seems to be related to the unsteady axial pinch, 
which is behind the origin of the Crab jet in this model. 
As we have mentioned already, the magnetic field in the 
backflow at the base of the jet is highly inhomogeneous and 
this results in strong spatial and temporal fluctuations of 
the magnetic hoop stress, and hence the axial pressure. As 
the result, the arch shock dives towards the pulsar at times 
when the pressure is high, and moves further out when it 
is low. Another factor is the presence of strong vortices in 
the backflow, which can appear all the way along the arch- 
shock. In their eyes the pressure is lower and on the outside 
it is higher. Moreover, there are significant fluctuations of 
the ram pressure inside the simulated nebula as well. 

In the simulations, the wisps are associated with regions 
of high magnetic field in the unsteady outflow from the ter- 
mination shock. New wisps are produced approximately on 
the light crossing time of the termination shock, which was 
around 10 months in the simulations and which is around 
3-4 months for the Crab Nebula. Strong variations of the 
shock structure occurred on the similar time scale. Thus, the 
variability of gamma-ray emission on the time-scale around 
several months may well have this origin. The variations of 
gamma-ray flux may be attributed to changes of the proper 
emissivity of the inner knot, associated with changes of the 
magnetic field strength and the number density of emitting 
particles, but also to changes in the direction of the Doppler 
beam (see Figure m). 

Assuming the power law spectral distribution for the 
emissivity in the comoving frame, j v oc u~ a , the observed 
emissivity is 



>,ob =D > 



( Lind fc Blandford|1985 1. When the angle between the flow 
direction and the line of sight decreases from <f> — I/7 to 
(f> — the Doppler factor increases from 7 to 27. Thus, the 
difference in the inclination angle of the arch shock at the 
location of inner knot A(j> — I/7 can bring about the dif- 
ference in the boosting factor up to 2 2+Q . Above 100 MeV, 
where the observed synchrotron spectrum can be approx- 
imated by a power law with a ~ 3 (Abdo et al. 20101, 



this corresponds to a 30-fold flux variation. At lower photon 
energies, where the knot spectral index is expected to be 
a ~ 0.6, the corresponding flux variation is five times lower. 
Moreover, the cooling time of the electrons, producing such 
photons, becomes significantly larger - it is already several 
years for the electrons emitting at 1 MeV. After traveling 
for such a long time the emitting plasma enters the remote 
parts of the nebula, where it inevitably decelerates and its 
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Figure 6. Variability of the gamma-ray emission from the inner knot in computer simulations jCamus et al.|2009fr . The five colour plots 
show the images of the inner knot at £ p h = 100 MeV (l v in linear scale). They are separated by 0.15 year, the time increasing from left 
to right and from top to bottom. The plot in the right bottom corner shows the corresponding total flux variation at £ p ^ = 100 MeV. 



emission is no longer subject to strong Doppler beaming 
The contribution of this unbeamed emission may explain 



the observations by de Jager et al. (19961, who noticed that 



during these observations the flux in the 1-30 MeV energy 
range was increasing, whereas in the 30-150 MeV range it 
was decreasing. Most likely, the total flux from the nebula 
at 1-30 MeV is dominated by few recently produced wisps. 

At 100 keV the observed total flu x from the Crab Neb- 
ula is ~ 4 x 10 _28 ergs _1 cm _2 Hz _1 (Horns & Aharonian 



|2004[ ), whereas the expected flux from the inner knot is 
-30„,„„-i„„-2xi„-i The lagt estimate is 

-0.64 



only ~ 1.3 x 10" du ergs _i cm"^Hz" 

obtained from the power law F„ oc u~ v o ' i , normalised using 
the observed optical flux of the inner knot ( |Tziamtzis et al.| 
20091. Thus, even a 10-fold increase of the X-ray emission 



from the knot would produce a variation below 5% in the 
total flux. For the similar reason, the knot variability could 
hardly be seen in 10 GeV-10 TeV range, which is dominated 
by the Inverse Compton emission of old electrons occupying 
the whole volume of the nebula ( de Jager et al.|1996 l. 

Finally, few words have to be said on the mysterious 
gamma-ray flares from the Crab Nebula. The observed linear 
sizes of the inner knot are l± — ip±d n ~ 6 light days and 
i|| = tp\\d n ~ 2 light days. This implies that the variability 



5 The observed deceleration of the Crab wisps is a clear confir- 
mation of such evolution I Hester et al.|2002 I . 



time-scale of the order of few days is possible, at least in 
principle. The fact that the flare spectrum extends beyond 
the highest characteristic energy allowed for the synchrotron 
emission can still be explained by the Doppler effect. Fermi 
reported the spectral index of the flaring component a ~ 1.5. 
Such a steep spectrum is expected because of the proximity 
of the exponential cut-off. The fact that flares are not seen 
at both lower and higher energies can be explained in the 
same fashion as for the long timescale variability. 

The most difficult task is to explain not only the short 
duration of the flares but also the fact that they are quite 
rare. If indeed they originated from the inner knot then this 
had to be associated with some rather peculiar events. For 
example, they could be produced when some explosive event 
at the base of Crab's jet drives a shock, which then collides 
with the termination shock near the knot location, with both 
shocks being almost parallel to each other. The shocks inter- 
section point could move towards the observer with speed 
very close to the speed of light, potentially resulting in a 
short burst of emission associated with this point. The mag- 
netic reconnection could be behind such explosions but this 
would probably require a significantly higher magnetization 
compared to what is assumed in the current MHD models. 
So far, the emission from simulated PWN was computed 
simply by integrating the instantaneous emissivity along the 
line of sight. The relativistic retardation effect was not taken 
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torus 



Figure 7. Variability of the termination shock as a reason behind 
the variability of high energy emission from the Crab Nebula. 
The solid and dotted lines show the termination shock with two 
extreme locations of its polar cusp. The dash-dotted lines show 
the corresponding directions of the Doppler beam. 



into account. This practice has to be abandoned in future 
studies as it filters out the short time-scale variability asso- 
ciated with the relativistic motion along the line of sight. 



5 CONCLUSIONS 

(i) Downstream of oblique termination shocks of pulsar 
winds, the Lorentz factor of bulk motion can be rather high 
even for low-sigma winds, up to 7 = 5 for reasonably small 
inclination angles. This is shown using the analytical solu- 
tion for oblique relativistic MHD shocks and confirmed by 
numerical simulations of the Crab Nebula. For a high-sigma 
wind the Lorentz factor can be even higher. 

(ii) The bright inner knot in the synthetic synchrotron 
maps of the Crab Nebula, obtained in relativistic MHD 
simulations, is definitely a highly Doppler-boosted emis- 
sion from the region located downstream of the termination 
shock (the so-called arch shock) and close to its polar cusp. 
The inner knot of the Crab Nebula is likely to be of the 
same origin. Its geometrical parameters, such as the ratio 
of its major axis to its distance from the pulsar, and the 
ratio of its minor and major axes, are consistent with this 
interpretation. 

(iii) The combination of the short synchrotron cool- 
ing time of gamma-ray emitting electrons and the strong 
Doppler beaming in the vicinity of the termination shock 
suggest that if the synchrotron gamma-ray electrons are ac- 
celerated mainly at the termination shock then the inner 
knot makes a major contribution to the integral gamma ray 
emission from the Nebula around 100 MeV. The two-point 
spectral index, a ~ 0.64, based on the observed optical emis- 
sion of the knot and the integral gamma-ray emission of the 
Crab Nebula at 100 MeV, is consistent with this hypothesis. 
A similar value is obtained in the "standard model" of the 



Crab Nebula emission by Kennel & Coroniti ( 1984 1 for the 



spectrum injected into the nebula by the termination shock. 

(iv) The observed variability of the Crab Nebula in the 
1-100 MeV window on the timescale from one month to sev- 
eral years can be related to the large scale variability of the 
termination shock discovered in recent high resolution nu- 
merical simulations. 



(v) The small size of the inner knot of the Crab Nebula, 2- 
6 light days, show that the recently discovered short gamma- 
ray flares can also originate from the knot. However, the 
exact mechanism behind such short and rare events remains 
unclear. 

The most critical prediction of our model, which allows 
a relatively simple test with currently available telescopes, is 
that the un-pulsed synchrotron gamma-ray emission of the 
Crab Nebula in the Fermi and AGILE windows originates 
from within one arcsecond of the pulsar itself. Although the 
angular resolution of gamma-ray telescopes is not even close 
to one arcsecond, the test can be based on comparing the 
gamma-ray light curve with the one obtained in optics for 
the inner knot. Within one arcsecond the inner knot is the 
dominant feature, apart from the pulsar itself. A potential 
problem of this test is the close proximity of the knot to the 
pulsar, which makes image based separation of their fluxes 
rather tricky, even for HST and the ground-based instru- 
ments with adaptive optics. However, this does not seem to 
be needed as the pulsed emission from the Crab pulsar is 
very stable and its un-pulsed emission is not expected to 
be variable too. Hence, one only needs to measure the total 
flux from within ~ 1" of the pulsar and subtract from it 
the phase averaged flux of the pulsed emission. The optical 
variability is expected to be strong, as the typical flux from 
the the inner knot is already about 6 — 10% of the flux from 
the pulsar flux. The light curves in optics and gamma rays 
should correlate. 
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APPENDIX A: RELATIVISTIC OBLIQUE MHD 
SHOCKS 

The first issue is how high can be the Lorentz factor down- 
stream of an oblique shock. In the limit of vanishing magne- 
tization, oblique relativistic shocks was studied by Konigl 
111980). The general case of relativistic MHD shocks was 



considered in Majorana & Anile (19871, whereas (Lyutikov 



20041 explored the special case where the magnetic field is 



parallel to the shock front. Here we deal only with this case, 
as the magnetic field of pulsar winds is almost perfectly asy- 
muthal and hence parallel to the axisymmetric wind termi- 
nation shock, and focus on the question of how relativistic 
can be the flow behind this shock. 

In the shock frame, the fluxes of energy, momentum, 
rest mass, and magnetic field are continuous across the shock 



(w + B 2 )f 2 /3 X = const, 
(w + B 2 )~f 2 l3 x (3 x +p + 



B 2 



const, 



(w + B )7 f3 x /3 y = const, 
P'fPx = const, 
B'yPx = const, 



(Al) 

(A2) 

(A3) 
(A4) 
(A5) 



where p is the rest mass density, p is the gas pressure, w = 
pc 2 + kP is the relativistic enthalpy, k = T/(T — 1), where V 
is the adiabatic index, B is the magnetic field as measured 
in the fluid frame, /3 = v/c, and 7 is the Lorentz factor. 
We select the frame where the velocity vector is in the xy- 
plane, the magnetic field is parallel to the z-direction, and 
the shock front is parallel to the yz-plane. In what follows 
we will use subscripts 1 and 2 to denote the upstream an 
the downstream states respectively. 



From Equations (All and (A3 1 it follows that 



Ply = I3 2 y, 



(A6) 



whereas Equations ( Al I , ( A4 1 and (A5l yield 



P2 _ £>2 _ CF2 
Pi til ""J 

1 

K 



V 

x' 



2 

a-2 



77(1 + Kdj + <7l) — <7l 



(A7) 
(A8) 



where a = B 2 /47r pc 2 is the magnetization parameter, a 2 
p/pc 2 is the temperature parameter, \ = fcx/fiix, and 



7l /1 1 2 /, 2nn1/2 
7)= — = (1 + W lx (l - X )) • 

72 



(A9) 



Equations ( A6)A9 1 allow to find the downstream state given 
the parameters of the upstream state and \. The equation 
for x is derived via substituting the expressions for p2, cr 2 , 
and a 2 from Equations ( A6| A8 \ into Equation (A2l, which 
can be written as 



(1 + khj +o"i) h pi(a ± + -ai) 

pi 2 



(1 + Hat. + G 2 )— + P2(«2 + ^2), 

P2 2 



(A10) 



where j = py/3 x - In general, this is a rather combersome al- 
gebraic equation which has to be solved numerically. How- 
ever, the pulsar winds are expected to be cold (ai — > 0) and 
highly relativistic (71 ^> 1), which allows significant simpli- 
fications. In the limit of cold upstream flow, this equation 
reduces to 



X 



(1 + o-i)u la 



-aiu\ x 



2 k 



+ 



(1 + ai)u lx 



1± 
2 



+ 



(1 + o-i)(1 + mL) 



A' 

v- 2 
2k 



(1 + u\ x ) 



+ 



= 0. 



(AH) 



Further simplification can be made if Ui x 2> 1. If we 
denote the angle between the shock plane and the velocity 
vector as <5, then this condition implies 81 2> l/7i- In this 
case, the above equation reduces to the cubic equation for x 

X S (1 + cti)(k-1)-x 2 (ai — 



X(l + ^i) 



+ K] + 







(A12) 



Its root x = 1 corresponds to the continious solution. Of the 
other two roots only 



X 



2 + 6a! + VV 



(A13) 



4(1 + err) (k- 1)' 
where 

V = (2 + aiK) 2 - 8(1 + ai)ffi(n - 1)(2 - k), 

is physical. For the ultrarelativistic equation of state, with 
r = 4/3 and k. = 4, this becomes 



X ; 



1 + 2<7i + ^16al + 16<ti + 1 
6(l + ai) ' 



(A14) 
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Figure Al. Downstream Lorentz factor of oblique shocks as a 
function of the shock angle. Solid lines correspond to the upstream 
Lorentz factor 71 = 10 6 and the magnetization parameter a = 
0.01, 1.0, and 10. The dashed line corresponds to 71 = 10 and 
cr = 0.01. 



The corresponding downstream Lorentz factor is 

72= , 1 ^V- (A15) 

For (7i ^> 1 this yields 

sin 81 
and for 

2\/2 2 'smc5i 



72 * 4h-> ( A16 ) 



Figure Al shows 72(^1) for o\ = 0, 1 and 10. 
The downstream shock angle 

tan $2 = x( <J i)tan5i. (A18) 

For a <C 1 and in the of small angle approximation this 
yields 

82 * ^1. (A19) 

The general conlusion from this analysis is that the 
downstream flow can be highly relativistic, provided the 
shock is sufficiently oblique. The current two-dimensional 
MHD models of PWN do predict a very oblique termina- 
tion shock due the anisotropic distribution of energy flux in 
the pulsar wind. For o\ <C 1 in the pulsar wind, the case 
favoured by the models, and in the small angle approxima- 
tion, we have 72 — l/<5i. The Lorentz factor can be signifi- 
cantly higher if a > 1. 

In the numerical simulations cited above the Lorentz 
factor of pulsar wind was rather low, 71 ~ 10. The results 
for such a relatively low upstream Lorentz factor are shown 
in Figure [Al] by the dashed line. One can see, that this really 
becomes a factor only for very small shock angles, 8\ < 0.2. 



